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With  the  funding  provided  by  AFOSP .  the  PI  was  able  to  take  his  sabbatical  lea%e 
for  the  period  of  1  Sep  91  -  31  Aug  92  at  Rutgers  University.  At  Rutgers,  the  FI 
collaborated  with  Prof.  Michael  Vogellus  on  several  problems.  The  research 
resulted  in  several  reports,  some  ot  which  have  already  been  accepted  for  ^ 

publication.  The  main  research  proposed  was  to  investigate  various  aspects  of  ^he 
problem  of  electrical  Impedan-  lomographv.  In  electrical  impedance  tomography, 
the  problem  is  to  image  the  interior  of  an  objecr  from  static  electrical 
measurements  conducted  on  the  boundary.  The  techniques  of  electrical  Impedance 
tomography  may  be  applied  to  nondestructive  evaluation  of  structures.  The  grant 
was  directly  responsible  for  two  roa ior  accomplishments  -  t 1)  Inversion  from  real 
data  and  (2)  Sensitivity  of  electrical  Impedance  tomography. 
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Summary 

With  the  funding  provided  by  AFOSK.  the  investigator  \i.as  able  to  take  his  sabbatif-ii  leave 
for  the  period  August  1991  to  .August  1992  at  Rutgers  fniversity.  .Ac  Rutgers,  tlie  pi  col¬ 
laborated  with  Prof.  Michael  Wtgelius  on  several  problenis.  The  research  resulted  in  se-.era.i 
reports,  some  of  which  have  already  been  accepted  for  publication. 

The  main  research  proposed  wa.s  to  investigate  various  aspert.s  of  the  problem  of  elerincaJ 
impedance  tomography.  In  electrical  impedance  tomography,  the  problem  is  to  image  the 
interior  of  an  object  from  static  electrical  measurements  conducted  on  the  boundary.  The 
techniques  of  electrical  impedance  tomography  may  be  applied  to  nondestructive  evaluation 
of  structures. 

The  grant  has  been  directly  responsible  for  two  major  arcornplisliments, 

1.  Inversion  from  real  data  In  collaboration  with  Prof.  \a!dis  I,i<’pa  of  th“  Radia¬ 
tion  Laboratory  at  the  University  of  .Michigan,  we  constructed  an  ekrtrical  impedance 
imaging  device  in  the  laboratory.  Data  collected  by  the  device  were  inverted  using  al¬ 
gorithms  that  Prof.  Vogelius  and  I  ha%e  developed  over  the  past  few  years.  The  restsits 
show  that  electrical  impedance  tomography  is  a  viable  alternative  to  weil-established 
nondestructive  evaluations  methods  such  as  eddy  current  and  ultrasonics, 

2.  Sensitivity  of  electrical  impedance  tomography  With  Prof.  David  Dobson  of  the 
School  of  Mathematics  at  University  of  .Minnesota,  we  developed  a  method  for  analyzing 
the  sensitivity  and  the  resolution  power  of  various  electrical  impedance  imaging  devices. 
The  results  further  leads  to  an  understanding  of  how  to  obtain  higher  resolution  images 
in  electrical  impedance  tomography. 

While  at  Rutgers,  the  PI  al.so  vvorked  with  Prof.  Abgelius  on  a  problem  arising  in  the 
study  of  vibrations  of  composite  materials.  We  developed  an  a.symptotic  theory  to  predict  the 
long-time  vibration  of  composite  with  a  periodic  microstnicture.  Our  re.suits  show  that  the 
long-time  behavior  of  the  vibration  depends  very  strongly  on  the  interaction  of  the  periodic 
structure  with  the  boundary  of  the  medium.  Rigorous  treatment  of  the  boundary  layer  effects 
on  the  vibration  was  carried  out.  The  work  resulted  in  two  reports,  a  third  is  being  written 
for  the  proceedings  of  a  conference  on  control  of  distributed  parameter  systems  whicli  was 
held  at  the  Institute  for  Mathematics  and  its  Applications  at  the  University  of  Minnesota. 
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Figure  1;  The  design  of  the  tank. 


Application  of  electrical  impedance  tomography  to  crack  de¬ 
tection 

We  set  out  with  the  goal  of  showing  that  the  method  of  electrical  impedance  tomography  is  a 
viable  alternative  to  the  problem  of  the  detection  (and  determination)  of  cracks  in  structures. 
To  this  end,  we  constructed  a  cylindrical  tank,  to  which  12  electrodes  are  attached.  Please 
consult  Figure  1.  The  tank  is  filled  with  an  electrolyte  solution.  The  presence  of  a  crack  is 
simulated  by  immersing  a  metal  strip  vertically  in  the  lank. 

To  generate  the  data,  a  low  frequency  .4C  source  is  used.  A  single  measurement  consists 
of  passing  current  through  an  adjacent  pair  of  electrodes  and  recording  the  resulting  voltage 
drops  in  all  adjacent  pairs  of  electrodes.  .A  total  of  12  measurements  is  carried  out.  cycling 
the  applied  current  through  all  possible  pairs  of  adjacent  electrodes. 

All  the  equipment  we  u.sed  —  current  source,  digital  voltmeter,  IEEE  buses,  etc;  are  of 
the  the  “off-the-shelf’  variety.  The  data  collected  were  then  processed  to  conform  with  the 
data  needed  for  the  two  algorithms  tested. 

The  two  algorithms  used  svere  designed  for  different  purposes.  The  first  one,  a  back- 
projection  algorithm,  was  first  proposed  by  Barber  and  Brown.  The  method  provides  an 
approximate  image  of  the  conductivity  distribution  in  the  interior  of  the  tank.  It  gives  only 
a  blurry  image  of  the  crack:  it  can  be  used  for  crack  detection  and  for  providing  a  rough 
estimate  of  the  crack  size  and  location. 

The  second  algorithm  has  been  specifically  designed  to  look  for  a  straightline  crack.  Given 
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Figure  2:  The  images  obtained  by  applying  the  back  project  ion  algorithm  on  the  rneavured 
data.  The  exact  locations  of  the  cracks  are  shown  for  reference. 


exact  data,  it  can  find  the  actual  location  of  the  crack  iteratively  beginning  with  an  initial 
guess. 

The  results  of  the  application  of  the  backprojection  algorithm  are  displayed  in  Figure.s  2a 
and  2b  for  two  experiment.  In  the  first  experiment,  the  actual  crack  is  a  line  segment  whose 
endpoints  are 

(-0.38,0.38)  and  (0.38.0.38). 

(in  units  of  the  cross  section  radius).  The  second  experiment  involve  a  2  inch  metal  strip 
simulating  a  crack  with  endpoints  at 

(0.24.0.50)  and  (0.62.0.50). 

In  each  figure,  we  superposed  ^he  actual  locations  of  the  crack  for  reference.  In  viewing  these 
results,  we  must  keep  in  mind  that  there  is  a  substantial  amount  of  noise  in  the  data  in 
the  form  of  measurement  errors,  as  well  a.s  some  modeling  errors.  These  results  do  however 
indicate  that  the  data  are  sufficiently  accurate  for  the  backprojection  algorithm  to  reconstruct 
a  rough  ims^e  of  the  cracks.  However,  the  images  are  too  blurry  to  give  good  estimates,  for 
instance  of  the  crack  sizes.  In  genercil,  even  with  perfect  data,  this  is  probably  the  best  we 
can  hope  for  with  the  backprojection  algortihm. 

The  result  of  application  of  the  crack  finding  algorithm  on  the  experimertal  data  is  given 
next.  In  applying  it  to  the  first  data  set.  we  took  the  initial  guess  to  be  a  crack  whose 
endpoints  are  (-0.5,  0)  and  (0.5.0).  The  algorithm  converged  in  4  iterations.  The  result  of 
this  calculation  is  shown  in  Figure  3.  which  gives  a  comparison  of  the  actual  crack  (shown 
in  dashes)  with  the  crack  found  by  the  algorithm.  The  reconstruction  yielded  a  crack  whose 
endpoints  are  (-0.3402,  0.3173)  and  (0.3885.  0.3341). 

An  initial  guess  in  the  form  of  a  crack  with  endpoints  at  (0.0)  and  (0.5.0)  was  used  when 
applying  the  algorithm  to  the  second  data  set.  Convergence  is  achieved  after  5  iterations.  A 
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Figure  3:  A  comparsion  of  the  actual  crack  with  the  crack  reonstructed  by  our  algorithm. 
The  algorithm  took  4  iterations  to  reconstruct  the  crack  shown  in  solid  line.  The  actual  crack 
is  shown  in  dashes. 


Figure  4:  A  comparison  of  the  actual  crack  with  the  crack  reconstructed  by  the  algorithm. 
The  algorithm  took  5  iterations  to  reconstruct  the  crack  shown  in  solid  line.  The  actual  crack 
is  shown  in  dashes. 
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comparison  of  the  actuaJ  location  of  the  crack  with  the  reconstructed  cr.o  k  i-  tn.sde  m  1  la  ;:. 
d,  The  reconstruction  yielded  a  crack  whose  end  points  are  •  (>  19t>-4,  d  }G»d>  <  assd  .  U.'.Tui. 
0.1773). 

Based  on  these  results,  we  propose  an  inspection  strategy  wliicii  cunsists  <jf  tlr't 
the  backprojection  algorithm  to  detect  the  presence  of  a  crack.  If  a  crack  is  detei  ted.  w,*  u-* 
the  result  from  the  backprojection  algorithm  to  obtain  an  initial  guess  f(jr  tlie  crack  fnidis-ir 
algorithm.  These  initial  findings  seem  to  indicate  that  FJectncal  Impedance  loruograph't 
may  prove  to  be  a  promising  nondestructive  evaluation  technnjue  for  certain  applicatioii' 
The  results  of  our  investigation  is  reported  in  publication  3  of  tlse  next  section,  which  ha"- 
been  submitted  for  publication  in  Journal  of  Nondestructive  Evaluations. 


Sensitivity  analysis  of  electrical  impedance  imaging  devices 

The  object  of  this  work  was  to  determine  that  a  certain  set  of  current  pat  tern  has  an  advaiitaije 
over  another  set  in  electrical  impedance  tomography.  .X  current  pattern  is  the  distribution 
of  the  applied  current  on  the  electrodes.  Barber  and  Brown  favor  the  “dipole"  pattern, 
and  Isaacson  and  his  co-worker,s  favor  the  “trigonometric"  pattern.  Besides  settling  th.e 
controversy  of  dipole  versus  trigonometric  patterns,  we  also  wanted  to  know  how  the  current 
pattern  is  responsible  for  an  increase  in  stability  and  resolution  in  an  electrical  impedance 
imaging  device. 

What  is  novel  in  our  work  is  that  we  do  not  discretize  the  problem.  Bather,  we  work 
with  a  concept  of  a  pseudo-inverse,  which  allow  us  to  conclude  properties  about  an  imaging 
device  without  discretization.  In  this  way.  we  were  able  to  conclude  convincingly  if  a  certain 
current  pattern  does  have  an  advantage  over  another.  The  quality  we  try  to  quantify  is  the 
stability  property  and  the  resolution  power  of  the  imaging  device. 

W^e  are  lead  to  the  conclusion  that  in  realistic  applications,  the  trigonometric  current 
pattern  is  better  than  the  dipole  pattern.  The  work  lead  to  a  direction  for  further  research: 
Current  patterns  that  focus.  If  we  are  interested  only  in  imaging  the  conductivity  of  a 
part  of  a  domain,  then  our  work  shows  that  it  is  possible  to  generate  a  .set  of  current  patterns 
that  gives  the  best  “illumination’'  of  that  part  of  the  domain.  This  is  a  subject  of  our  current 
research.  The  work  described  here  has  appeared  as  an  I.MA  (Institute  for  Mathematics  and 
its  Applications)  report,  and  has  been  submitted  to  SI.\M  Journal  of  Applied  Mathematics 
for  publication. 
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1.  First  order  corrections  to  the  homogenized  eigenvalues  of  a  periodic  composite  medium 
(with  M.  Vogelius).  to  appear  in  SIAM  J.  Appl.  Math. 
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3.  Crack  determination  from  boundary  measurements  —  reconstruction  from  experimental 
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